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Abstract 

The  computer  language  Logo  facilitates  the  teaching  of  analytic  geometry  and  calculus  from  the  notion  of 
curvature,  through  its  "turtle  geometry"  facility  {2].  We  provide  some  theoretical  basis  for  finding  turtle 
geometry  equivalents  of  familiar  curves  in  analytic  geometry,  and  vice  versa,  by  some  simple  methods 
apparently  previously  unnoticed.  In  particular,  we  study  turtle  geometry  programs  where  the  curvature  of  a 
line  is  a  trigonometric  function  of  its  orientation. 

1.  Introduction 

Recently  the  computer  language  Logo  has  become  available  for  the  first  time  on  a  wide  variety  of  personal 
computers  [1, 6, 3].  A  distinctive  feature  of  Logo  is  an  emphasis  in  the  graphics  package  on  "turtle  geometry", 
a  geometry  emphasizing  discrete  changes  to  local  curvature  of  curves  rather  than  coordinate  systems.  Despite 
its  simplicity,  turtle  geometry  is  quite  powerful,  and  provides  an  alternative  viewpoint  that  makes  certain 
mathematical  concepts  clearer,  as  discussed  in  [2], 

But  an  issue  not  well  addressed  in  [2]  is  relationship  of  curves  created  by  turtle  geometry  to  the  familiar 
ones  of  analytic  geometry.  In  this  paper,  thus,  we  present  methods  for  finding  small-step  turtle-geometry 
equivalents  (mostly  previously  unnoticed)  for  some  well-known  analytic-geometry  curves,  and  vice  versa. 
This  material  provides  a  good  exercise  for  calculus  students,  while  also  demonstrating  some  numerical 
analysis  concepts. 

2.  Turtle  geometry 

Turtle  geometry  draws  planar  curves  using  two  commands,  FORWARD  and  RIGHT.  A  cursor  called  a 
turtle  represents  a  point  and  orientation  on  a  graphics  screen.  FORWARD  causes  the  turtle  to  move  forward 
in  a  straight  line  a  distance  (approximately  in  millimeter  units)  specified  as  its  argument,  in  the  direction  the 
turtle  is  pointing.  RIGHT  causes  the  turtle  to  turn  in  place  a  number  of  degrees  specified  by  its  argument 

Turtles  thus  have  a  state  consisting  of  a  screen  position  and  an  orientation  (or  heading).  The  Cartesian 
coordinates  are  referred  to  as  XCOR  and  YCOR,  and  the  heading  as  HF.AD1NG  (though  we  will  use  the 
abbreviation  H  for  it  in  this  paper).  Wc  assume  the  turtle  starts  with  XCOR=0.  YCOR=0  (the  center  of  the 
screen),  and  HEADING=0  (pointing  straight  up).  And  a  HEADING  of  90  is  straight  to  the  right,  180 
straight  down,  and  270  straight  to  the  left  Note  that  HEADING  can  be  thought  of  as  the  total  amount  a 
turtle  has  turned  right  since  it  started. 

Logo  programs  arc  written  in  terms  of  procedures.  A  procedure  consists  of  a  sequence  of  commands,  on 
successive  lines,  'flic  first  line  contains  the  word  TO  followed  by  the  name  of  the  procedure,  and  then 
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arguments  if  any.  'Hie  last  line  of  a  procedure  is  signalled  by  END. 

3.  Curvature 

We  start  r-om  the  classical  definition  of  curvature  in  a  Cartesian  coordinate  system,  as  given  in  most 
advanced  calculus  treatments  (e.g.  [4],  p.  305),  assuming  that  clockwise  is  positive  curvature: 

K(x,y)  =  -(d2y/dx2)/[l  +  (dy/dx)2]15 

We  can  put  this  in  a  form  closer  to  turtle  geometry  by  noting  that  dy/dx  =  cot(H),  H  the  heading  of  the  turtle 

[5J: 

K(x,y,H)  =  -(d2y/dx2)  /  [1  +  cot2H]15  =  -(d2y/dx2)  *  sinJH 
We  call  this  the  "hybrid"  curvature  formula  because  it  includes  both  Cartesian  and  heading  variables. 

Curvature  is  the  rate  at  which  a  curve  turns  clockwise,  and  thus  can  be  approximated  in  Ix>go  by  the  ratio 

of  the  amount  a  turtle  turns  RIGHT  to  the  amount  it  goes  FORWARD  --  provided  the  turtle  takes  small 

enough  steps  so  as  to  make  its  path  look  sufficiently  smooth.  (Practically  speaking,  this  means  something  like 

steps  of  1  millimeter  and  turtle  turns  of  five  degrees  on  most  graphics  terminals  that  use  Logo.)  So  if  we 

define  the  curvature  function  we  want  a  curve  to  follow  by  the  procedure  K  of  no  arguments,  cither  of  these 

recursive  Logo  procedures  will  probably  draw  what  we  want: 

TO  CURVE 
FORWARD  1 
RIGHT  K 
CURVE 
END 

TO  CURVE 

FORWARD  (QUOTIENT  l  K) 

RIGHT  1 

CURVE 

END 

We  shall  call  them  the  "angle-control"  and  "step-control"  variants,  respectively.  Which  should  we  use? 
Unfortunately,  sometimes  they  do  not  work  the  same,  and  only  one  gives  to  the  correct  result  This  is  because 
these  arc  discrete  approximations  of  smooth  curves,  and  anytime  the  absolute  value  of  K  becomes  large  in  the 
first  version,  or  the  absolute  value  of  K  approaches  0  in  the  second,  the  approximation  breaks  down.  So  we 
will  have  to  watch  for  that. 

In  what  follows  we  shall  only  concern  ourselves  with  the  shape  of  a  curve,  and  not  its  size  or  orientation. 
Thus  we  shall  ignore  constant  multipliers  in  front  of  a  curvature  formula,  except  sometimes  the  sign,  and  we 
shall  throw  away  any  constant  multipliers  generated  in  analysis  routinely.  Wc  also  will  not  concern  ourselves 
with  at  what  screen  location  wc  start  drawing  the  shape  (except  in  section  4.3),  since  turtle  geometry  docs  not 
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use  a  coordinate  system,  though  the  usual  convention  in  Logo  is  to  start  at  the  center  of  a  graphics  screen.  But 
it  will  make  a  difference  what  heading  we  start  with  in  the  procedures  we  discuss  --  different  pieces  of  die 
curve  will  appear  for  different  starting  headings,  and  for  some  headings  not  defined  for  the  curve  (c.g. 
asymptotes),  pathological  curves  like  straight  lines  will  be  generated. 

Special  note:  in  this  paper  we  assume  that  headings  of  greater  than  360  degrees  or  less  than  0  are 
meaningful  and  significant  as  indications  of  "how  far  the  turtle  has  turned".  Most  Logo  implementations, 
however,  take  "heading"  as  this  modulo  360.  As  a  result,  some  of  the  Logo  procedures  discussed  below  will 
not  execute  properly  in  those  implementations  (in  particular,  those  where  heading  (H)  appears  other  than  as 
an  argument  to  a  trigonometric  function).  But  it  is  easy  to  define  a  counter  variable  to  use  instead  of 
HEADING. 

4.  Some  examples:  from  y  =  f(x)  to  K  =  g(H) 

We  give  now  some  examples  of  curves  initially  described  as  y  =  f(x)  in  new  terms  of  their  curvature  as 
function  of  heading,  K  =  g(H).  The  method  is  (1)  find  dy/dx  and  dfydx2  in  terms  of  x;  (2)  express  d*y/dx2 
in  terms  of  dy/dx  only,  no  x’s;  and  (3)  then  substitute  cot(H)  for  dy/dx,  and  multiply  by  -sinJ(H)  to  get  the 
curvature  K  in  terms  of  H.  This  method  is  not  guaranteed  to  work  because  step  (2)  is  not  possible  for  some 
fix),  but  it  will  work  for  some  interesting  functions. 

4.1 .  Conic  sections 

Let  us  apply  the  hybrid  formula  to  some  familiar  curves.  First,  the  general  parabola  y  =  ax2  +  bx  + 

c.  Since  dfydx2  =  2a,.a  constant,  following  our  policy  of  ignoring  constants  we  just  turn  right  sinJ(H)  at  each 

turn  in  the  "angle-control"  procedure,  as  follows: 

TO  PARABOLA 
FORWARD  1 

RIGHT  (PRODUCT  (SIN  HEADING)  (SIN  HEADING)  (SIN  HEADING)) 

PARABOLA 

F.ND 

If  the  turtle  is  started  with  a  heading  between  0  and  180,  it  will  draw  a  downward-tending  parabola;  else,  an 
upwards-tending  one. 

One  of  the  first  things  people  learn  to  do  with  Logo  is  to  draw  a  circle  by  repeating  many  times  a 

forward-right  pair  of  commands  with  constant  arguments: 

TO  CIRCLE 
FORWARD  1 
RIGHT  1 
CIRCLE 
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END 

(This  is  a  special  case  of  the  well-known  Logo  procedure,  POLY.)  We  can  prove  that  in  fact  y  =  -/(l^-x2)  is 
drawn  with  this  procedure,  dy/dx  =  -x/V( r^x2),  and  d2y/dx2  =  rVf^-x2)15.  Hence  1  +  (dy/dx)2  = 
r2/(r2-x2),  and(l  +  (dy/dx)2)1 5/r  =  d2y/dx2.  Substituting  cot(H)  =  dy/dx,  d2y/dx2  =  (1  +  cot2(H))15  = 
csc3(H).  Hence  the  curvature  is  -sin3(H)  *  csc3(H)  =  -1,  a  constant 

We  can  generalize  this  to  ellipses  and  hyperbolas.  First  consider  the  general  ellipse  y  =  c -/(a2-x2),  where  c 

is  the  ratio  of  the  height  to  the  the  width.  Then  dy/dx  =  -cx/-/(a2-x2)  and  d2y/dx2  =  ca2/(a2-x2)is.  Hence 

1  +  (dy/dx)2/c2  =  a2/(a2-x2),  and  c(l  +  (dy/dx)2/c2)15  =  d2y/dx2  =  c(l  +  (l/c2)cot2(H))15  =  (l/c2Xc2 

+  cot2(H))15.  Multiplying  by  -sin3(H),  and  ignoring  the  multiplicative  constant  this  then  leads  to  a  total 

curvature  of  -(c2sin2(H)  +  cos2(H))3  5  =  -[(c2-l)sin2(H)  +  l]15.  Since  sin2(H)>0  and  c2-l>-l  for  c*0,  the 

curvature  is  always  negative  for  c*0,  but  the  degree  of  curve  varies  periodically  with  the  heading.  Here  is  the 

Logo  procedure  for  this  general  ellipse: 

TO  ELLIPSE  :C 
FORWARD  1 

RIGHT  THRF.EHALVES  DIFFERENCE  (PRODUCT  (DIFFERENCE  SQUARE  :C  1) 

SQUARE  SIN  HEADING)  1 

ELLIPSE  :C 
END 

TO  THREEHALVES  :X 
SQRT  (PRODUCT  :X  :X  :X) 

END 

By  similar  analysis  we  can  get  the  curvature  of  a  general  hyperbola,  y  =  cV(a2+x2).  It  turns  out  dy/dx  = 
cx/V(a2+x2)  and  d2y/dx2  =  ca2/(a2+x2)15,  so  c(l  -  (l/c2)cot2(H))’  5  =  d2y/dx2,  and  the  curvature  is 
proportional  to  -(c2sin2(H)  -  cos2(H))15  =  *{(c2+l)sin2(H)  -  1J15.  This  has  maximum  curvature  at  H=90 
degrees,  and  then  gradually  turns  left  more  and  more  slowly,  reaching  an  asymptote  when  the  curvature 
becomes  zero,  or  (c2+l)sin2(H)  =  1,  orH  =  arcsin[/(l/(l+c2)]  =  arccot(c). 

Note  both  the  ellipse  and  hyperbola  become  a  parabola  when  c=0,  since  (l-sin2(H))!  5  =  cos3(H).  And 
when  c  approaches  oo,  both  approach  (c2sin2(H))’ 5  =  c3sin3(H),  which  is  also  a  parabola. 

4.2. Inverses 

Take  l(x)=e\  dy/dx  =  d2y/dx2  =  e\  so  d2y/dx2  =  cot(H),  and  the  curvature  needed  to  draw  it  is 
-sin3(H)  *  cot(H)  =  -sin2(H)cos(H)  =  *.5sin(l!)sin(2H). 

Since  the  inverse  function  of  fix)  looks  just  like  it  but  "turned  over”  and  oriented  in  a  different  direction. 


we  should  expect  their  curvature  formulas  to  be  related.  Wc  can  verify  this  for  the  inverse  function  of  e\ 
ln(x);  dy/dx  =  1/x  and  d2y/dx2  =  -1/x2,  so  d2y/dx2  =  -cot2^,  and  the  curvature  is  sin(H)cos2(H)  =  .S 
cos(H)sin(2H). 

In  general,  we  can  obtain  the  curvature  of  the  inverse  of  a  function  by  interchanging  all  sin(H)  and  cos(H) 
terms  and  then  reversing  the  sign  of  the  total  expression  so  that  it  curves  in  the  opposite  direction.  This  is 
because 

dx/dy  =  1  /  (dy/dx)  =  tan(H) 

So  tan(H)  should  be  substituted  for  cot(H)  in  all  occurrences  of  it,  and  vice  versa.  But  since 

1  +  (dy/dx)2  =  1  +  cot2(H)  =  l/sin2(H) 

1  +  (dx/dy)2  =  1  +  tan2(H)  =  l/cos^H) 

sin(H)  and  cos(H)  must  necessarily  interchange  too. 

4.3.  Problems  with  inflection  points 

Unfortunately,  wc  usually  cannot  use  an  expression  for  the  curvature  in  terms  of  only  the  heading  when  the 
curvature  changes  sign  somewhere.  Often  these  conditions  arise  with  square  roots  which  can  be  cither 
positive  or  negative.  An  example  is  y  =  x3;  dy/dx  =  3x2,  and  d2y/dx2  =  6x.  So  trying  to  express  d2y/dx2  in 
terms  of  dy/dx,  we  get  ±  V(12(dy/dx)]  =  ±  V[12cot(H)J.  and  a  curvature  of  =Psm3(H)  V[12cot(H)]. 

But  we  cannot  choose  only  one  sign  to  cover  the  whole  range  --  the  sign'  must  change  at  the  inflection  point, 

the  point  where  6x  =  0,  x=0,  or  the  place  where  arccot[3(0)2]  =  H  =  90.  We  must  therefore  add  a 

conditional  branch,  something  like  this  (assuming  wc  start  with  a  heading  which  is  a  small  positive  number): 

TO  CUBE 
FORWARD  1 

RIGHT  PRODUCT  (SQUARE  SIN  HEADING)  (SQRT  COT  HEADING) 

IF  HEADING  =  90  [CUBERIGHTSIDE]  ELSE  [CUBE] 

END 

TO  CUBERIGHTSIDE 
FORWARD  I 

RIGHT  MINUS  PRODUCT  (SQUARE  SIN  HEADING)  (SQRT  COT  HEADING) 

CUBERIGHTSIDE 

END 

where  SQRT  is  the  square  root  function  (see  the  Logo  manuals  cited  for  the  exact  interpretation  of  this 

syntax).  But  since  the  curve  is  really  composed  of  small  straight  lines,  with  some  roundoff  error  in  calculating 

the  heading  over  a  period  of  time,  wc  may  never  actually  reach  a  heading  of  7,cro.  Instead  wc  had  better  say 

for  procedure  CUBE: 

TO  CUBE 
FORWARD  1 
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RIGHT  PRODUCT  (SQU ARF.  SIN  HEADING)  (SQRT  COT  HEADING) 

IF  HEADING  <  89  (CUBE)  ELSE (CUBERIGHTSIDEJ 

END 

with  CUBERIGHTSIDE  remaining  as  above. 


5.  Parametric  equations:  from  x  =  f^t),  y  s  f2(t)  to  K  =  g(H) 

Some  curves  are  easier  to  express  in  terms  of  parametric  equations  for  x  and  y  than  in  terms  of  a  equation 
including  both  x  and  y.  To  find  curvature  in  terms  of  heading  for  them  there  are  two  steps:  (1)  express  die 
parametric  equations  in  terms  of  H  entirely  (use  dy/dx  =  cot(H)),  and  (2)  compute  derivatives  with  respect  to 
this  expression  to  obtain  the  curvature.  We  can  rewrite  the  basic  curvature  formula  in  terms  of  derivatives  of 
H: 

K(H)  =  -sin3(H)*d2y/dx2 

But  d2y/dx2  =  d/dx  (dy/dx)  =  d/dx  (cot(H))  =  csc2(H)  *  dH/dx 
Hence  K(H)  =  -sin(H)  /  dx/dH 

But  also 

dH/dx  =  dH/dy  *  dy/dx  =  dH/dy  *  cot(H) 
so  also  K(H)  =  -cos(H)  /  dy/dH 


As  an  example,  take  the  cycloid  described  by 
x  =  a(t  -  sin(t)),  y  =  a(l  -  cos(t)) 

cot(H)  =  dy/dx  =  (dy/dH)/(dx/dH)  =  sin(t)  /  (1  -  cos(t))  =  cot(t/2) 

and  hence  t/2  =  H,  t  =  2H,  and 

x  =  a(2H  -  sin(2H)),  y  =  a(l  -  cos(2H)) 
dx/dH  =  2a(l  -cos(2H)),  dy/dH  =  2asin(2H) 

So  the  curvature  is,  using  the  simpler  expression  dy/dH: 

-sin(F'  '  2a  sin(2H)  =  -1  /  4a  cos(H) 

and  hence  the  figure  can  be  drawn  by  the  Logo  procedure 

TO  CYCLOID 

FORWARD  SIN  HEADING 

RIGHT  1 

CYCLOID 

END 

Note  we  use  the  "step-control"  procedure  format  because  scc(H)  becomes  infinite  at  ±  180  degrees. 


6.  Polar  coordinates  to  K  s  g(H) 

We  can  also  draw  curves  given  in  polar  coordinates,  r  =  f(ff): 
x  =  rcos(0),  y  =  rsin(0) 

dx/d0  =  (dr/d0)cos(0)  -  r  sin(0),  dy/d0  =  (dr/dfl)sm(0)  +  rcos(fl) 
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cot(H)  =  dy/dx  =  (dy/dff)/(dx/d$)  = 

[(dr/dtf)sin(0)  +  r  cos(0)]  /  [(dr/d0)cos(0)  -  r  sin(0)j 
=  cos(tf  -  Arctan[(dr/d0)/r]]  /  -sin[0  -  Arctan[(dr/d0)/r]] 
=  -cot[0  -  Arccoit.,  (dr/d#)]] 

Hence 

H  =  -6  +  Arccot(r/(dr/d#)],  or  cot(H  +  6)  =  r/(dr/dfl) 


The  curvature  of  a  curve  expressed  in  polar  coordinates  can  be  found  by  manipulation  of  the  x-y  curvature 

formula,  and  is  given  in  many  books  (e.g.  [4j,  p.  306): 

K(6)  =  [r2  +  2(dr/d#)2  -  r(d2r/d#2)]  /  [r2  +  (dr/dfl)2)ls 

We  substitute  into  it  the  relationship  just  derived,  i.e.  (dr/d#)/r  =  tan(H+0): 

K(#,H)  =  [1  +  2tan2(H + 0)  -  (d2r/d#2)/r]  /  ifl  +  tan2(H+#)]ls 
=  [2cos(H  +  0)  -  cos3(H + 0)  -  [cos3(H  +  6)  *  (d2r/d#2)/r]]  /  r 
=  [2cos(H  +  0)\/x  -  [cos3(H +  #)*[!  +  (dVdflfyr2]) 


For  particular  situations  we  may  be  able  to  express  the  above  only  in  terms  of  H.  For  example,  take  the 
logarithmic  spiral  r  =  ccB.  for  which  tan(H+0)  =  (dr/d#)/r  =  ce^/e6*  =  c.  Hence: 

dr/d#  =  ccc*  d2r/d# 2  =  c2ec* 

K(H,r)  =  2cos(arctan(c))/r  -  cos3(arctan(c))  *  fl  +  c2J/r 

Which  just  factors  into  a  constant  term  divided  by  r  =  e0*  =  ec<arccot<c>'I0  Hence  the  curvature  is  just  a 

constant  times  e'cH  =  b‘H,  b = ec.  So  it  can  be  drawn  by  this  procedure: 

TO  I.OGSPIRAL 

FORWARD  EXPONENT  1.01  H 

RIGHT  l 

LOGSPIRAL 

END 

where  EXPONENT  represent  the  first  argument  taken  to  the  power  represented  by  the  second  argument. 


7.  Changing-argument  recursive-procedure  forms 

There's  a  simpler  way  we  can  express  I.OGSPIRAL  (given  in  chapter  2  of  [2]),  by  giving  it  an  argument: 

TO  LOGSPIRAL  :SIDE 
FORWARD  :SIDE 
RIGHT  1 

LOGSPIRAL  (:SIDE*  1.01) 

END 

where  :SIDE  is  a  parameter  to  the  procedure. 


Many  procedures  that  draw  smooth  curves  can  be  written  in  this  form.  Letting  the  curvature  on  two 
successive  recursive  calls  be  Kj  and  K2,  we  try  to  find  a  simple  expression  relating  K2  to  K,,  preferably 


independent  of  heading.  In  the  above.  K2 = K  j/1.01  which  is  independent  of  heading.  Note  if  the  curve  is 
smooth  enough,  we  can  approximate  K2  by  a  Taylor  series  about  Kj,  that  is 
K,(H)  =  KjIH)  +  (5K(H)/5H)AH 

so  there  is  always  an  "additive"  approximation  for  any  Logo  step-control  procedure: 

TO  CURVE  .CURVATURE 
FORWARD  1 
RIGHT  CURVATURE 

CURVE  SUM  CURVATURF.  (DF.RIVK  HEADING) 

END 

where  DF.RIVK  is  a  procedure  that  computes  the  derivative  of  the  curvature  function.  But  since  this  is  only 
an  approximation  it  will  (usually)  gradually  diverge  from  the  true  curve,  though  slower  for  some  functions 
than  others. 


8.  From  K  =  g(H)  to  y  =  f(x) 

Now  we  address  the  inverse  problem,  which  is  to  determine  the  Cartesian  description  of  a  function  given 
its  curvature  as  a  function  of  heading.  Two  approaches  are  possible. 

8.1 .  From  g(H)  to  y  =  f(x) 

l  et  /.  =  dy/dx.  Then  if  the  curvature  is  a  function  of  heading  g(H): 

g(H)  = -sin5(H)*dz/dx 
-csc3(arccol(z))  *  g(arccot(/))  =  dz/dx 
dx  =  d/  /  (  csc'(arccol(/))  *  g(arecol(z))] 
x  =  -Jd//(g(arccot(7)Xl  +  z2)15l 

In  particular  eases  we  may  be  able  to  solve  this,  then  invert  it  and  integrate  with  respect  to  x  to  get  an  fl[x). 

For  instance,  consider  the  curve  drawn  by  the  procedure 

TO  MYSTFRY1 
FORWARD  1 
RIGID  SIN  HEADING 
MYSTERY  1 
END 

which  draws  a  kind  of  "squashed  parabola".  Its  curvature  is  K(H)  =  sin(t  l).  Hence  the  equation  is 

x  =  -/d//(l  F  /?)  -  C  -  arctan(z) 

So  /.  -  un(C-x)  =  dy/dx. 

Hence  y  =  /  un(C-x)dx  =  ln|cos(C-x)| 

which  is  thus  the  mystery  curve  .as  a  function  of  x. 

As  another  example,  consider  the  similar  procedure 
TO  MYSTKRY2 
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FORWARD  1 

RIGHT  PRODUCT  SIN  HEADING  SIN  HEADING 

MYSTERY2 

END 

which  draws  a  less-squashed  parabola.  Its  curvature  is  K(H)  =  sin2(H),  and: 

x  =  -/dz/-/([l+z2]  =  C-arcsinh(z) 

So  z  =  sinh(C-x)  =  dy/dx. 

Hence  y  =  /  sinh(C-x)dx  =  cosh(C-x) 

which  is  a  catenary. 


This  method  depends  however  on  three  steps,  two  integrals  and  a  function  inversion.  If  any  of  these  are 

not  possible,  the  method  will  not  work.  But  anytime  the  first  integral  is  possible,  one  docs  derive  a  differential 

equation  in  terms  of  dy/dx  and  x,  and  it  may  be  possible  to  solve  this  with  other  methods,  or  perhaps 

approximate  a  solution  by  infinite  series.  For  instance,  the  procedure 

TO  MYSTERY3 
FORWARD  1 

RIGHT  EXPONENT  (SIN  HEADING)  5 

MYSTERY3 

END 

leads  to  the  differential  equation 
(dy/dx)3/3  +  (dy/dx)  =  x 
which  has  no  easy  solution. 


8.2.  From  g(H)  to  x  =  f(H)  and  y  =  h(H) 

If,  however,  it  is  acceptable  to  determine  the  curve  in  terms  of  the  parameter  H,  then  we  can  often  solve 
other  situations  too.  Recall  the  formulas  from  section  5  for  curvature  in  terms  of  the  heading  parameter 
exclusively: 

K(H)  =  -sin(H)  /  (dx/dH) ,  or  K(H)  =  -cos(H)  /  (dy/dH) 

We  can  rearrange  these  expressions  for  K(H)  and  integrate  them: 

dx  =  -sin(H)dH  /  K(H).  dy  =  -cos(H)dH  /  K(H) 
x  =  -/[sin(H)dlI  /  K(H)] 
y  =  -/  (cos(H)dll  /  K(H)] 

Sometimes  we  can  then  find  a  substitution  to  eliminate  H  and  get  an  equation  in  terms  of  only  x  and  y,  but 
not  always. 

As  an  example,  consider  curves  drawn  by  the  "step-control"  form 
TO  AMPOLY  :C 

FORWARD  SIN  PRODUCT  :C  HEADING 
RIGHT  1 


♦ 
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AMPOLY  :C 
END 

where  :C  denotes  a  parameter.  (AMPOLY  stands  for  "amplitude  modulation  POLY",  since  the  curvature  of 

a  circle  is  being  "modulated"  at  some  frequency.)  Here  the  curvature  is  l/sin(cH),  and  hence 

x  =  Jsin(H)sin(cH)dH  =  ,5/{cos((l-c)H)  -  cos((l+c)H)}dH 
=  (l/2(l-c))sin((l-c)H)  -  (1/2(1 +c))sin((l+c)H) 
y  =  /cos(H)sin(cH)dH  =  .5J[sin((l+c)H)-sin((l-c)H)]dH 
=  (l/2(l-c))cos((l-c)H)  -  (1/2(1 +c))cos((l  +c)H) 

These  are  parametric  equations  for  general  hypocycloids  and  epicycloids  ([4],  p.  312-314),  with  the 
substitutions  t/2  =  H  and  (l-c)/2  =  b/a).  When  c>l,  it  is  a  hypocydoid;  when  c<l,  an  epicycloid;  and  when 
c= 1,  a  cycloid,  as  we  showed  in  section  S.  We  can  ignore  the  constants  of  integration  --  they  just  affect  the 
location  of  the  shape. 


As  another  example,  take  the  Logo  procedure  ( [2],  ch.  1) 

TO  POLYSPI  :S1DE 
FORWARD  :SIDE 
RIGHT  1 

POLYSPI  SUM  :SIDE  1 
END 

which  can  also  be  written,  following  the  discussion  of  section  7,  and  assuming  a  starting  heading  of  0,  as 

TO  POLYSPI 
FORWARD  HEADING 
RIGHT  1 
POLYSPI 
END 

Hence  K(H)  =  -1/H.  Then  from  integration  by  parts: 

x  —  /  Hsin(II)dH  =  sin(H)  -  Hcos(H) 
y  =  /  Hcos(H)dH  =  cos(H)  -  Hsin(H) 

We  cannot  find  an  easy  nonparamctric  form  in  general,  but  note  as  H  approaches  infinity  the  second  terms 
predominate,  so 

x  25  -Hcos(H),  y  25  -Hsin(H) 
r  =  V[x2  +  y2]  =  H,  6  =  arctan[tan(H)]  =  H 

and  r=0  is  a  polar  expression  for  the  curve.  Thus  as  H  becomes  large  the  curve  approaches  the  spiral  of 
Archimedes  ( [4],  p.  317).  Wc  can  show  in  general  that 
TO  SPIRAL  :N 

FORWARD  EXPONENT  HEADING  :N 
RIGHT  1 
SPIRAL  :N 
END 

approaches  r=0N  as  H  becomes  large,  for  NX). 


8.3.  Connections  to  Fourier  analysis 

The  inverse  of  the  curvature  K  is  often  called  the  "radius  of  curvature"  R.  Our  equations  for  parameterized 

x  and  y  can  be  written  as 

x  =  -/  R(H)sin(H)dH 
y  =  -/  R(H)cos(H)dH 

and  one  can  see  these  are  just  the  formulas  for  the  coefficients  of  sin(H)  and  cos(H)  in  a  Fourier  series 
expansion  of  R(H). 

There  is  another  connection  to  Fourier  analysis,  however.  Note  if  R(H)  =  R^H)  +  R2{H),  then 

x  =  -/  R,(H)sin(H)dH  -  /  R2(H)sin(H)dH 
y  =  -/  Rj(H)cos(H)dH  -  /  R2(H)cos(H)dH 

and  the  resulting  curve  is  the  "vector  addition"  of  the  two  original  curves,  as  in  DUOPOLY  of  ch.  3  of  [2J. 
Now  since  we  can  approximate  any  curve  sufficiently  accurately  with  a  sufficient  number  of  terms  from  its 
Fourier  series,  we  can  approximate  an  arbitrary  R(H)  as  the  sum  of  weighted  sines  and  cosines: 

R(H)  =  ZJacosOH/c)  +  b^infiH/c)! 

where  c  is  a  constant  representing  how  fine  accuracy  we  desire  in  the  approximation.  Hence  we  can  express  x 
and  y  as 

x  =  -2/(acos(iH/c)+bsin(iH/c)lsin(H)dH 
y  =  -2i/[a.cos(iH/c)+bjsin(iH/c)lcos(H)dH 

and  the  result  is  a  weighted  vector  sum  of  hypocycloids  and  epicycloids,  in  the  same  way  that  MULTIPOLY 
in  ch.  3  of  [2]  is  the  sum  of  circles  of  different  sizes  and  curvatures. 

9.  Conclusion 

We  have  shown  a  variety  of  methods  for  obtaining  small-step  turtle  geometry  programs  from  equations  for 
analytic-geometry  curves,  and  vice  versa.  Unfortunately,  each  method  only  works  part  of  the  time,  and  no 
guarantees  can  be  given.  Sec  [5]  for  further  student  problems  related  to  this  material. 
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